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1 Introduction 



One-dimensional harmonic oscillator H^ H °^ = p 2 + x 2 is exactly solvable and appears 
in virtually any textbook on quantum mechanics. In the recently proposed VT 
symmetric quantum mechanics of Bender et al |I| a similar guiding role is played by 
the non-Hermitian cubic Hamiltonian H^ co ^ = p 2 +i x 3 . A straightforward numerical 
and semi-classical analysis of its Schrodinger equation 

H^\^ n )=E^\^ n ), n = 0,l,... (1) 

supports the highly plausible conjecture that the spectrum of energies is real, discrete 
and bounded below 0. The conjectured absence of its imaginary components is also 
indicated by the Hilbert-Schmidt analysis || and by the perturbation calculations 
in both the weak-coupling regime [[| and its strong-coupling, purely numerically 
generated re-arrangements ||. Also the numerical analytic continuation |J and the 
so called exact semiclassical analysis J!]] reveal numerous parallels to the Hermitian 
models. 

Definitely, the problem deserves a deeper study. In what follows we shall modify 
the traditional perspective and work in the momentum representation. It gives the 
momentum as a mere number, p G (—00,00), and represents the coordinate x by 
the operator x = id p . Equation ([!]) acquires a new form containing the purely real 
differential Schrodinger operator H^ co ^ — E of the third order. This gives an unusual 
formulation of our bound-state problem, 

^(p) = E^(p) , ^(p) e L 2 (M). 

The quadratic p— dependence of the kinetic term T(p) = p 2 does not seem to make 
the equation any easier to solve. Here, we shall modify the kinetic term and offer 
the comprehensive analysis of the resulting square-well-type approximate form of the 
imaginary cubic oscillator. 
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2 Local solutions in a solvable simplification 



In a way proposed by Priifer || many wave functions can be visualized as certain 
"deformations" ip(x) ~ cisin[^(x)] + C2Cos[g(x)] of solutions which correspond to a 
locally constant potential in the standard quantum mechanics. Such a type of an 
ansatz finds immediate applications in numerical computations || and leads to an 
easy interpretation of the traditional Sturm Liouvillean oscillation theorems [11|. 
Similar idea will also be used here as our main methodical guide. 



2.1 Square- well approximation of the kinetic term 

Let us replace the kinetic energy operator T(p) = p 2 by the most elementary square 
well of a finite depth Z > 0, 



T{p) 



Z, p 6 (-oo,-l), 

o, 

Z, pe(l,oo). 



In the bounded range of energies E < Z this enables us to split our toy model in the 
two separate differential equations, 
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iP(p) = 0, 



P e (-oo,-l) |J(l,oo) 



The two auxiliary parameters a = a(E) > and f3 = (3(E) > are defined in 
such a way that Z = E + 8/3 3 >E = 8 a 3 , and appear in the three independent 
(exponential) solutions of our equations. Their general superpositions are complex 
but may also be given the real, trigonometric form. In this way, near the origin we 
shall write 



%{p) = de 2ap + fe- ap cos(ap + 



pe (-1,1), 
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where the symbols d, f and 9 stand for the three undetermined constant parameters. 
In the right and left asymptotic regions we obtain the similar formulae. After we omit 
their unphysical and normalization- violating, exponentially growing components, we 
get the one-parametric family to the right, 

Mp) = ge~ Wp , pe(i,oo). 

Its two-parametric left-barrier counterpart reads as follows, 

ii;-(p) = ce? p cos(j3p + r)), p<E(-oo,-l), p = V3/3 . 

At the right discontinuity p = 1 we have to guarantee the continuity of ip(p), d p ip(p) 
and dpi/)(p). This is equivalent to the three matching conditions 

de 2a + f e~ a cos(a + 9) = g e~ 2/3 , 
2ad e 2a -af e~ a [cos(a + 9) + Vs sin(a + 9)] = -2(3 g e" 2/3 , (2) 
Aa 2 d e 2a - 2a 2 f e~ a [cos(a + 9)-^ sin(a + 9)} = 4(3 2 g e~ 213 . 
The weighted sum of these equations re-scales and interrelates the unknown coeffi- 
cients d and g in terms of a new energy parameter t = t(E) = (3(E) /a(E) > or 
rather R = R(E) = (1 - t + t 2 y x l 2 > 0, 

de2a = D{E) = ^wy G w=9< r * p - 

After we eliminate G from the last two equations (0) which are linear in d and / we 
obtain the elementary formula which defines the shift 9 = 9(E), 

tanfa + 9) = — ; . 

v ' 2/t-l 

The trigonometric factors become fixed up to their common sign e — ±1, 

/3 

cos(a + 9) = e [1 - t(E)/2] R(E), sm(a + 9) = ^— e t(E) R(E). 
The same sign also enters our last decoupled definition 
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Thus, up to an overall normalization (say, g = 1), all the free parameters in our wave 
functions matched at p = 1 become specified as functions of the energy. 

At the second discontinuity p = — 1 we have to satisfy the other three matching 
conditions. In terms of the abbreviations 

L x = d e~ 2a + fe a cos(-a + 0), 
L 2 = 2ade~ 2a - a f e a [cos(-a + 6) + V% sin(-a + 0)], 
L 3 = Aa 2 d e~ 2a - 2a 2 f e a [cos(-a + 9) - V% sm(-a + 9)] 

they read 

Li(a, d, f, 9) = ce"* 3 cos(— (3 + rj), 
L 2 (a,d, f,9) = (3 c e~ P {cos(- f3 + ri) - sin(-/3 + 77), 
L 3 (a, d, f, 9) = -2p 2 ce- (3 [cos(-f3 + 77) + \/3sin(-/? + 77). 

In the same manner as above they determine the values of c = c(E) and 77 = rj(E). 
Moreover, their properly weighted sum gives the equation 

1+1+2 

y/3F(E) sin(-a + 9) + (2i - 1) cos(-a + 9) + 2 ^— .0(E) e~ 6a = 

or, after a slight re- arrangement, the amazingly transparent relation 

(l-4t + t 2 ) cos(2 v / 3a) + v / 3(l-t 2 ) sin^v^a) = e ~ 3 °) ■ ( 3 ) 

It should be read as our desired final implicit definition of the physical energies E. 

2.2 Z— dependence of the toy spectrum 

It is quite instructive to search for the physical energies numerically. Starting from 
the very-shallow-well extreme in eq. we find the two clearly distinguished energy 
roots. The qualitative features of the graph remain unchanged in a broad interval 
of the "inverse strengths" 1/Z. Its shape is sampled in Figure 1 at Z = 1CT 3 . We 
have tested that even the approximate height w —5 of its left plateau stays virtually 
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unchanged in between Z = and Z = 10~ 3 . Within the same interval of the 
shallowest wells the left zero grows from the value 0.00280 till 0.0241. Beyond the 
broad, downwards-oriented peak one finds the second, right zero moving from the 
value 0.01047 (found very close to the instantaneous threshold 0.01077) up the value 
0.1056 (not very far from its threshold 0.1077, either) within the same interval of 
1/Z. 

A new feature emerges around Z = 10 _1 (with the left zero at 0.0445 and with 
the right zero 0.222 still quite close to the threshold 0.232) and Z = 1 (with the left 
zero at 0.072 and with the right zero 0.446, not that close already to the threshold 
1/2) in the left half of the picture. The plateau develops a local, safely negative 
maximum. 

In the subsequent domain of Z > 1 we have to switch our attention back to the 
right half of our graph. Immediately before the coupling reaches the integer value 
Z — 5, the end of the curve returns to the negative half-plane near the maximal (i.e., 
threshold) energy. This means that there emerges the third energy level there. The 
total number of bound states grows to N = 3 (cf. the leftmost items in our Table 1). 

Beyond Z — 5, our attention has to return quickly to the left half of the picture 
where the very slow growth of the local maximum creates a new quality at last. The 
top of the local bump touches and crosses the horizontal axis at Z pa 5.3003 and 
E pa 0.6244. At Z — 5.3005 a new doublet of energies is formed in a way illustrated 
in Figure 2. The number of levels jumps to N = 5. 

A smooth deformation of the graph takes place when the value of Z grows on. 
During this evolution we discover that our (originally broad), downward-oriented 
peak shrinks quite quickly and moves comparatively slowly to the right. It gets close 
to the rightmost and, to its bad luck, slightly more slowly moving zero number five. 
The magnified picture of the resulting "collision" is displayed here in Figure 3. At 
Z = 35 it shows that in the threshold region of the energies, 
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• the wavy motion of the threshold end of our graph still did not manage to 
reach the zero axis; 

• the downwards-oriented peak has already left the positive part (and moved to 
the negative part) of the curve in question. 

As a consequence, the number of levels drops, quite unexpectedly, down to 3 again 
(cf. Table 1). 

In the vicinity of Z = 40 the new, rightmost energy root emerges at last. Up 
to Z = 100 and beyond, the number of levels stays equal to 4. Then it increases 
to 5, due to the emergence of the next threshold zero. Only after that, the slowly 
moving downward peak reaches the domain of the fourth zero. At almost exactly 
Z = 190 its left (and temporarily negative) local maximum reaches the zero value 
again (cf. Figure 4). At this moment the number of states jumps up by two to seven. 
The magnified graphical proof is offered by Figure 5 at Z = 200. The latter Figure 
illustrates nicely the rapid shrinking of our peak with Z. The numerical detection of 
its position becomes more and more difficult. Although this position plays a crucial 
role in the practical determination of the number of levels at a given Z, we must 
be very careful in distinguishing the subgraph of Figure 5 (with three zeros) from a 
simple straight line with the single zero. 

The pattern is deceitful and the standard software which searches for roots has 
to be used with a due care. Vice versa, the above analysis enables us to take into 
account all the specific features of the Z— dependence of the graph in eq. (|j). We 
get a regular pattern summarized in Table 1 and exhibiting a certain regularity of 
the Z— dependence of the number of levels N = N(Z). 

3 Discussion 
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3.1 Energies in the square-well approximation 

The main "physical" consequence of the presence of the above-mentioned narrow 
peak is an unusual irregularity observed in the emergence of the new levels in the 
deeper wells. We may conclude that this irregularity is not an artifact of the com- 
putation method. The energy formula (|3|) for our square-well toy model is exact 
and the seemingly unpredictable emergence of its roots just reflects the fact that our 
Schrodinger equation is of the third order. In particular, there exists no symme- 
try/antisymmetry with respect to the parity p — > — p etc. 

Methodical consequences of our analysis are a bit discouraging. Firstly, the very 
symbolic-manipulation derivation of our present formulae proved unexpectedly com- 
plicated even in comparison with the multiple standard square wells in textbooks. 
That's why we did not move to any further piece- wise constant approximations 
of T{p). 

Secondly, even our use of the most elementary solvable example revealed quite 
clearly a very real danger of the possible loss of certain levels. For an illustration 
let us imagine that our numerical study would have been started in the deep-well 
domain, i.e., at the large Z. It is quite easy to generate the graphs of eq. (§) there. 
In the standard and routine finite-precision computer arithmetics one discovers that 
the results are very smooth and look virtually the same, say, in the interval of 
Z E (1000,1200). 

Let us pick up, for definitness, the larger sample Z = 1200. We get a picture 
(cf. our last Figure 6) which is regular and, deceptively, indicates that iV(1200) = 7. 
Unfortunately, the correct answer (appearing also at the right end of our Table 1) is 
iV(1200) = 10. Its derivation requires the use of a significantly enhanced precision. 
Otherwise, whenever we use just the standard 14 digits and Figure 6, we would have 
missed as much as three (i.e., cca 30 % of all) energy levels. 

In the light of our preceding considerations, an easy explanation of the latter 
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numerical paradox lies in the presence of the narrow peak. A priori, it is hardly 
predictable of course. It is necessary to spot it by brute force. One finds that at 
Z = 1000, this anomalous peak still lives safely below the sixth energy level. The 
related number of levels is reliably confirmed as equal to seven, indeed. 

In between Z = 1000 and Z = 1200, it is necessary to work in an enhanced 
precision arithmetics. One finds that the upper, threshold end of the curve crosses 
the horizontal axis only slightly above Z = 1100. Due to the very steep slope of the 
curve in this region, this crossing is not visible even at z — 1200 in Figure 6. 

One also has to trace the narrow peak carefully. It overtakes the sixth energy 
level at Z 1190 (and E = 4.217), in an arrangement resembling our Figure 4 
above. Thus, one concludes, finally, that the new, almost degenerate pair of the 
energy levels emerges immediately beyond this point. 

3.2 Wave functions and their zeros 

A marginal merit of our use of the square- well-shaped T(p) lies in the availability of 
the explicit wave functions. For the lack of space we have to omit all the illustrative 
pictures and mention just a few of their most characteristic features. 

In the first step we can notice that in the rightmost interval of p the absence of any 
nodal zero in the wave function is in fact very similar to the usual Sturm Liouvillean 
behaviour. Less expectedly, at the exact energy value one encounters an infinity of 
the nodal zeros in the leftmost subinterval of p G (— oo, — 1). In this domain we are 
fortunate in studying the exactly solvable case. The very presence of this infinite 
"left" set of nodal zeros is extremely sensitive to the numerical level of precision we 
use. Indeed, the errors are proportional to the unphysical ^( un P h y s -) (p) ~ exp(— 2(3 p) 
which is growing rapidly at p <C — 1 . 

We may also note that after the smallest deviation of the energy E from its abso- 
lutely precise bound-state value even the non-numerical and absolutely precise wave 
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functions will be dominated by the growing asymptotics x(j( un P h y s -)(j)^ ~ exp{— 2/3 p) 
near the left infinity of course. After all, the change of sign of these asymptotics 
remains to be a reliable source of information about the fact that the energy just 
crossed it physical value. In practice, this observation could survive in the appro- 
priate modifications of the "shooting" numerical algorithms [[TIJ or in the rigorously 
proven versions of the so called method of Hill determinants [|l"3|j). 



4 Outlook 

We may summarize that in the standard momentum representation, we represented 
our Hamiltonian with the imaginary cubic interaction by the real differential expres- 
sion. On a suitable Hilbert space this specifies the Hamiltonian operator with the 
numerical range (and, hence, spectrum) which is, obviously, real. This complements 
the extensive discussion of this topic in the recent literature (cf., e.g., || for its 
updated list). 

We have touched here several immediate constructive consequences of the latter 
important observation. Being interested in the detailed structure of the energy spec- 
trum first of all, we arrived at a consistent picture of some of its les standard features 
by using the suitable approximation of the kinetic term T(p). 

Our main observation is the regularity of the Z— dependence of the number N(Z) 
of the bound states (cf. Table 1). This indicates that one cannot rely upon any type 



of a simply modified Sturm-Liouvillean theory at present [13]. A slightly better 
chance of a possible future deduction of some simple oscillation-type theorems or 
rules seems only to exist within the middle interval of p G (—1,1). There, for a 
continuously growing energy parameter E, a steady right-ward movement of the 
nodal zeros competes with the exponential terms which are slowly varying. One can 
expect a further progress achieved by the further intensive study of the similar VT 



symmetric examples [ITH . 
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Our numerical experiments also inspire a number of the further open questions. 
Some of them emerge in the purely numerical context. In particular, a better under- 
standing of the behaviour of the nodal zeros in the wave functions could help us to 
formulate an appropriate generalization of the Priifer-type algorithms. One would 
like to find some rules replacing the standard oscillation theorems. Especially in the 
vicinity of the correct physical energies they could really lead, say, to some reliable 
and robust "right-to-left shooting" numerical recipe etc. 

Several qualitative aspects of the problem became clarified even within our rough 
piece-wise constant approximation of T(p). We saw how the behaviour of the wave 
function asymptotics differs in the left and right infinity. In the vicinity of the 
origin p « the emergence and motion of the nodal zeros was interpreted in a 
graphical manner explaining some features of the N(Z) dependence. The loss of its 
monotonicity was confirmed by our solvable model. 

The use of the momentum representation proved able to throw a new light on 
the counterintuitive bound states in VT symmetric quantum mechanics. The emer- 
gence/disappearance of our quasi-degenerate doublets should be emphasized as, per- 
haps, analogous to the unavoided level crossings in harmonic oscillators |TB[ and/or 



to the anomalous doubling of levels in the models of Natanzon type 17 ] . Perhaps, 
all the similar irregularities in the spectra should be attributed to some common 
peculiar combination of the analyticity and non-Hermiticity in all these unusual VT 
symmetric systems. 
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Tables 



Table 1. 



Number of levels N and its changes A with growing Z. 



N 


235 34578 679 10 


A 


12-21121-2121 
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Figure captions 
Figure 1. 

Graphical solution of equation @ at Z — 1/1000. 
Figure 2. 

Local maximum giving the new doublet of roots at Z — 5.3005. 
Figure 3. 

The local maximum not giving the doublet of roots at Z — 35. 
Figure 4. 

Typical Z > 1 graph of eq. (|) (Z = 190). 
Figure 5. 

Quasi- degeneracy of the doublet of roots at Z = 200. 
Figure 6. 

Numerical invisibility of the narrow peak and of the new threshold root at Z = 1200. 
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Figure 3. 
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1 Introduction 



One-dimensional harmonic oscillator H^ H °^ = p 2 + x 2 is exactly solvable and appears 
in virtually any textbook on quantum mechanics. In the recently proposed VT 
symmetric quantum mechanics of Bender et al [1] a similar guiding role is played by 
the non-Hermitian cubic Hamiltonian if( co ) = p 2 -\-ix 3 . A straightforward numerical 
and semi-classical analysis of its Schrodinger equation 

HW\ip n ) = E^\ip n ), n = 0,l,..- (1) 

supports the highly plausible conjecture that the spectrum of energies is real, discrete 
and bounded below [2]. The conjectured absence of its imaginary components is also 
indicated by the Hilbert-Schmidt analysis [3] and by the perturbation calculations 
in both the weak-coupling regime [4] and its strong-coupling, purely numerically 
generated re-arrangements [5]. Also the numerical analytic continuation [6] and the 
so called exact semiclassical analysis [7] reveal numerous parallels to the Hermitian 
models. 

Definitely, the problem deserves a deeper study. In what follows we shall modify 
the traditional perspective and work in the momentum representation. It gives the 
momentum as a mere number, p G (—00,00), and represents the coordinate x by 
the operator x = id p . Equation (1) acquires a new form containing the purely real 
differential Schrodinger operator H( co ^> — E of the third order. This gives an unusual 
formulation of our bound-state problem, 

iP(p) = EiP(p) , v(p) e L 2 (M). 

The quadratic p— dependence of the kinetic term T(p) = p 2 does not seem to make 
the equation any easier to solve. Here, we shall modify the kinetic term and offer 
the comprehensive analysis of the resulting square-well-type approximate form of the 
imaginary cubic oscillator. 
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2 Local solutions in a solvable simplification 

In a way proposed by Priifer [8] many wave functions can be visualized as certain 
"deformations" ip(x) ~ c\ sin[^(x)] + c 2 cos[^(a;)] of solutions which correspond to a 
locally constant potential in the standard quantum mechanics. Such a type of an 
ansatz finds immediate applications in numerical computations [9] and leads to an 
easy interpretation [10] of the traditional Sturm Liouvillean oscillation theorems [11]. 
Similar idea will also be used here as our main methodical guide. 

2.1 Square- well approximation of the kinetic term 

Let us replace the kinetic energy operator T(p) = p 2 by the most elementary square 
well of a finite depth Z > 0, 

z, pe (-00,-1), 
t(p) = \ o, pe(-i,i), 
z, pe(i,oo). 

In the bounded range of energies E < Z this enables us to split our toy model in the 
two separate differential equations, 

d 3 
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Hp) = 0, 



iP(p) = 0, 



pe (-1,1), 



pe (-00,-1) U(i,oo). 



The two auxiliary parameters a = a(E) > and j3 = f3(E) > are defined in 
such a way that Z = E + 8f3 3 >E = 8 a 3 , and appear in the three independent 
(exponential) solutions of our equations. Their general superpositions are complex 
but may also be given the real, trigonometric form. In this way, near the origin we 
shall write 



ip (p) = de 2ap + fe- ap cos(ap + 9), pe (-1,1), 
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where the symbols d, f and 9 stand for the three undetermined constant parameters. 
In the right and left asymptotic regions we obtain the similar formulae. After we omit 
their unphysical and normalization- violating, exponentially growing components, we 
get the one-parametric family to the right, 

i; + (p) = ge- 2 ^, pe (l,oo). 

Its two-parametric left-barrier counterpart reads as follows, 

t/)_(p) = ce Pp cos((3p + r)), pe(-oo,-l), /3 = \/3/5. 

At the right discontinuity p = 1 we have to guarantee the continuity of tp(p), d p tp(p) 
and d 2 tp(p). This is equivalent to the three matching conditions 

de 2a + f e- a cos(a + 9) = ge^, 
2ad e 2a -af e~ a [cos(« + 9) + \/3 sin(a + 9)] = -2(3 g e" 2/3 , (2) 
Aa 2 de 2a - 2a 2 fe~ a [cos(« + 6) - sin(« + 6)] = 4/?V~ 2/3 - 
The weighted sum of these equations re-scales and interrelates the unknown coeffi- 
cients d and g in terms of a new energy parameter t = t(E) = f3(E)/a(E) > or 
rather R = R(E) = (1 - t + t 2 )" 1 / 2 > 0, 

de 2 * = D(E) = ^y y G(E)=ge- 2 ?. 

After we eliminate G from the last two equations (2) which are linear in d and / we 
obtain the elementary formula which defines the shift 9 = 9(E), 

a/3 

tan(a + 9) = — ; . 

V ' 2/t-l 

The trigonometric factors become fixed up to their common sign e = ±1, 

cos(« + 9) = e [1 - t(E)/2] R(E), sin(« + 9) = — s t(E) R(E). 

2 

The same sign also enters our last decoupled definition 
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Figure 1 . 







Figure 1: Graphical solution of equation (3) at Z = 1/1000. 

Thus, up to an overall normalization (say, g = 1), all the free parameters in our wave 
functions matched at p = 1 become specified as functions of the energy. 

At the second discontinuity p = — 1 we have to satisfy the other three matching 
conditions. In terms of the abbreviations 

L x = de~ 2a + / e a cos(-a + 9), 
L 2 = 2ade~ 2a - a f e a [cos(-a + 0) + \/3 sin(-a + 0)], 
L 3 = 4a 2 d e- 2a - 2a 2 f e a [cos(-a + 9) - sin(-a + 9)] 

they read 

Li(a, d, f, 9) = ce _/3 cos(—f3 + 77), 
L 2 (a, d, /, 9) = /?ce _/3 [cos(-/5 + 77) - \/3sin(-/5 + 77), 
L 3 (a, d, f, 9) = -2/3 2 ce- /3 [cos(-/3 + 77) + \/3sin(-/3 + 77). 
In the same manner as above they determine the values of c = c(.E) and rj = Tj(E). 
Moreover, their properly weighted sum gives the equation 

1 - t + t 2 

V3F(E) sin(-Q + 9) + (2t - 1) F(E) cos(-a + 0) + 2 ^— £>(£) e~ 6a = 

or, after a slight re- arrangement, the amazingly transparent relation 

(l-4t + t 2 ) cos(2v / 3«) + v / 3(l-t 2 ) sm(2V3a)= ( ^"^V 3 ") • (3) 



4 



Figure 2: Local maximum giving the new doublet of roots at Z = 5.3005. 

It should be read as our desired final implicit definition of the physical energies E. 
2.2 Z— dependence of the toy spectrum 

It is quite instructive to search for the physical energies numerically. Starting from 
the very-shallow-well extreme in eq. (3) we find the two clearly distinguished energy 
roots. The qualitative features of the graph remain unchanged in a broad interval of 
the "inverse strengths" \jZ. Its shape is sampled in Figure 1 at Z = 10~ 3 . 

We have tested that even the approximate height m — 5 of its left plateau stays 
virtually unchanged in between Z = 10~ 5 and Z = 10~ 3 . Within the same interval 
of the shallowest wells the left zero grows from the value 0.00280 till 0.0241. Beyond 
the broad, downwards-oriented peak one finds the second, right zero moving from 
the value 0.01047 (found very close to the instantaneous threshold 0.01077) up the 
value 0.1056 (not very far from its threshold 0.1077, either) within the same interval 
of l/Z. 

A new feature emerges around Z = 10 _1 (with the left zero at 0.0445 and with 
the right zero 0.222 still quite close to the threshold 0.232) and Z = 1 (with the left 
zero at 0.072 and with the right zero 0.446, not that close already to the threshold 
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Figure 3: The local maximum not giving the doublet of roots at Z = 35. 



1/2) in the left half of the picture. The plateau develops a local, safely negative 
maximum. 

In the subsequent domain of Z > 1 we have to switch our attention back to the 
right half of our graph. Immediately before the coupling reaches the integer value 
Z = 5, the end of the curve returns to the negative half-plane near the maximal (i.e., 
threshold) energy. This means that there emerges the third energy level there. The 
total number of bound states grows to N = 3 (cf. the leftmost items in our Table 1). 

Beyond Z = 5, our attention has to return quickly to the left half of the picture 
where the very slow growth of the local maximum creates a new quality at last. The 
top of the local bump touches and crosses the horizontal axis at Z m 5.3003 and 
E ?x 0.6244. At Z = 5.3005 a new doublet of energies is formed in a way illustrated 
in Figure 2. The number of levels jumps to N = 5. 

A smooth deformation of the graph takes place when the value of Z grows on. 
During this evolution we discover that our (originally broad), downward-oriented 
peak shrinks quite quickly and moves comparatively slowly to the right. It gets close 
to the rightmost and, to its bad luck, slightly more slowly moving zero number five. 
The magnified picture of the resulting "collision" is displayed here in Figure 3. At 
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Figure 4 . 









Figure 4: Typical Z > 1 graph of eq. (3) (Z = 190). 

Z = 35 it shows that in the threshold region of the energies, 

• the wavy motion of the threshold end of our graph still did not manage to 
reach the zero axis; 

• the downwards-oriented peak has already left the positive part (and moved to 
the negative part) of the curve in question. 

As a consequence, the number of levels drops, quite unexpectedly, down to 3 again 
(cf. Table 1). 

In the vicinity of Z = 40 the new, rightmost energy root emerges at last. Up to 
Z = 100 and beyond, the number of levels stays equal to 4. Then it increases to 5, 
due to the emergence of the next threshold zero. Only after that, the slowly moving 
downward peak reaches the domain of the fourth zero. At almost exactly Z = 190 
its left (and temporarily negative) local maximum reaches the zero value again (cf. 
Figure 4). 

At this moment the number of states jumps up by two to seven. The magnified 
graphical proof is offered by Figure 5 at Z = 200. 
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Figure 5: Quasi-degeneracy of the doublet of roots at Z = 200. 

The latter Figure illustrates nicely the rapid shrinking of our peak with Z. The 
numerical detection of its position becomes more and more difficult. Although this 
position plays a crucial role in the practical determination of the number of levels at 
a given Z, we must be very careful in distinguishing the subgraph of Figure 5 (with 
three zeros) from a simple straight line with the single zero. 

The pattern is deceitful and the standard software which searches for roots has 
to be used with a due care. Vice versa, the above analysis enables us to take into 
account all the specific features of the Z— dependence of the graph in eq. (3). We 
get a regular pattern summarized in Table 1 and exhibiting a certain regularity of 
the Z— dependence of the number of levels N = N(Z). 
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Figure 6: Numerical invisibility of the narrow peak and of the new threshold root at 
Z = 1200. 

3 Discussion 

3.1 Energies in the square- well approximation 

The main "physical" consequence of the presence of the above-mentioned narrow 
peak is an unusual irregularity observed in the emergence of the new levels in the 
deeper wells. We may conclude that this irregularity is not an artifact of the com- 
putation method. The energy formula (3) for our square-well toy model is exact 
and the seemingly unpredictable emergence of its roots just reflects the fact that our 
Schrodinger equation is of the third order. In particular, there exists no symme- 
try/antisymmetry with respect to the parity p —> —p etc. 

Methodical consequences of our analysis are a bit discouraging. Firstly, the very 
symbolic-manipulation derivation of our present formulae proved unexpectedly com- 
plicated even in comparison with the multiple standard square wells in textbooks. 
That's why we did not move to any further piece- wise constant approximations 
of T(p). 
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Secondly, even our use of the most elementary solvable example revealed quite 
clearly a very real danger of the possible loss of certain levels. For an illustration 
let us imagine that our numerical study would have been started in the deep-well 
domain, i.e., at the large Z . It is quite easy to generate the graphs of eq. (3) there. 
In the standard and routine finite-precision computer arithmetics one discovers that 
the results are very smooth and look virtually the same, say, in the interval of 
Z E (1000,1200). 

Let us pick up, for definitness, the larger sample Z = 1200. We get a picture 
(cf. Figure 6) which is regular and, deceptively, indicates that iV(1200) = 7. Un- 
fortunately, the correct answer (appearing also at the right end of our Table 1) is 
iV(1200) = 10. Its derivation requires the use of a significantly enhanced precision. 
Otherwise, whenever we use just the standard 14 digits and Figure 6, we would have 
missed as much as three (i.e., cca 30 % of all) energy levels. 

In the light of our preceding considerations, an easy explanation of the latter 
numerical paradox lies in the presence of the narrow peak. A priori, it is hardly 
predictable of course. It is necessary to spot it by brute force. One finds that at 
Z = 1000, this anomalous peak still lives safely below the sixth energy level. The 
related number of levels is reliably confirmed as equal to seven, indeed. 

In between Z = 1000 and Z = 1200, it is necessary to work in an enhanced 
precision arithmetics. One finds that the upper, threshold end of the curve crosses 
the horizontal axis only slightly above Z = 1100. Due to the very steep slope of the 
curve in this region, this crossing is not visible even at z = 1200 in Figure 6. 

One also has to trace the narrow peak carefully. It overtakes the sixth energy 
level at Z m 1190 (and E = 4.217), in an arrangement resembling our Figure 4 
above. Thus, one concludes, finally, that the new, almost degenerate pair of the 
energy levels emerges immediately beyond this point. 
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3.2 Wave functions and their zeros 

A marginal merit of our use of the square- well-shaped T(p) lies in the availability of 
the explicit wave functions. For the lack of space we have to omit all the illustrative 
pictures and mention just a few of their most characteristic features. 

In the first step we can notice that in the rightmost interval of p the absence of any 
nodal zero in the wave function is in fact very similar to the usual Sturm Liouvillean 
behaviour. Less expectedly, at the exact energy value one encounters an infinity of 
the nodal zeros in the leftmost subinterval of p e (— oo, — 1). In this domain we are 
fortunate in studying the exactly solvable case. The very presence of this infinite 
"left" set of nodal zeros is extremely sensitive to the numerical level of precision we 
use. Indeed, the errors are proportional to the unphysical ip( un P h y s -)(p} ~ exp(— 2(3 p) 
which is growing rapidly at p <C — 1. 

We may also note that after the smallest deviation of the energy E from its abso- 
lutely precise bound-state value even the non-numerical and absolutely precise wave 
functions will be dominated by the growing asymptotics ip( un P h y s -) (p) ^ exp(— 2(3 p) 
near the left infinity of course. After all, the change of sign of these asymptotics 
remains to be a reliable source of information about the fact that the energy just 
crossed it physical value. In practice, this observation could survive in the appro- 
priate modifications of the "shooting" numerical algorithms [12] or in the rigorously 
proven versions of the so called method of Hill determinants [13]). 

4 Outlook 

We may summarize that in the standard momentum representation, we represented 
our Hamiltonian with the imaginary cubic interaction by the real differential expres- 
sion. On a suitable Hilbert space this specifies the Hamiltonian operator with the 
numerical range (and, hence, spectrum) which is, obviously, real. This complements 
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the extensive discussion of this topic in the recent literature (cf., e.g., [3] for its 
updated list). 

We have touched here several immediate constructive consequences of the latter 
important observation. Being interested in the detailed structure of the energy spec- 
trum first of all, we arrived at a consistent picture of some of its les standard features 
by using the suitable approximation of the kinetic term T(p). 

Our main observation is the regularity of the Z— dependence of the number N(Z) 
of the bound states (cf. Table 1). This indicates that one cannot rely upon any type 
of a simply modified Sturm-Liouvillean theory at present [14]. A slightly better 
chance of a possible future deduction of some simple oscillation-type theorems or 
rules seems only to exist within the middle interval of p € (—1,1). There, for a 
continuously growing energy parameter E, a steady right-ward movement of the 
nodal zeros competes with the exponential terms which are slowly varying. One can 
expect a further progress achieved by the further intensive study of the similar VT 
symmetric examples [15]. 

Our numerical experiments also inspire a number of the further open questions. 
Some of them emerge in the purely numerical context. In particular, a better under- 
standing of the behaviour of the nodal zeros in the wave functions could help us to 
formulate an appropriate generalization of the Priifer-type algorithms. One would 
like to find some rules replacing the standard oscillation theorems. Especially in the 
vicinity of the correct physical energies they could really lead, say, to some reliable 
and robust "right-to-left shooting" numerical recipe etc. 

Several qualitative aspects of the problem became clarified even within our rough 
piece- wise constant approximation of T(p). We saw how the behaviour of the wave 
function asymptotics differs in the left and right infinity. In the vicinity of the 
origin the emergence and motion of the nodal zeros was interpreted in a 
graphical manner explaining some features of the N(Z) dependence. The loss of its 
monotonicity was confirmed by our solvable model. 
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The use of the momentum representation proved able to throw a new light on 
the counterintuitive bound states in VT symmetric quantum mechanics. The emer- 
gence/disappearance of our quasi-degenerate doublets should be emphasized as, per- 
haps, analogous to the unavoided level crossings in harmonic oscillators [16] and/or 
to the anomalous doubling of levels in the models of Natanzon type [17]. Perhaps, 
all the similar irregularities in the spectra should be attributed to some common 
peculiar combination of the analyticity and non-Hermiticity in all these unusual VT 
symmetric systems. 
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Number of levels N and its changes A with growing Z. 



N 


235 34578 679 10 


A 
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